Measurement of scintillation efficiency for nuclear recoils in liquid argon 
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The scintillation light yield of liquid argon from nuclear recoils relative to electronic recoils has 
been measured as a function of recoil energy from 10 keVr up to 250 keVr at zero electric field. 
The scintillation efficiency, defined as the ratio of the nuclear recoil scintillation response to the 
electronic recoil response, is 0.25 ± 0.02 + 0.01 (correlated) above 20 keVr. 



PACS numbers: 

I. INTRODUCTION 

A number of existing and proposed experiments use 
liquefied noble gases as detection media for Weakly In- 
teracting Massive Particles (WIMPs) [11-01 > ^ moti- 
vated dark matter candidate Q. Liquefied noble gases 
have a high scintillation yield, are relatively simple to 
purify of both radioactive contaminants and light ab- 
sorbers, and should be easily scalable to the large masses 
required for very sensitive detectors. Although the best 
limit for the spin-independent WIMP-nucleon cross sec- 
tion is currently set by the germanium-based CDMS ex- 
periment at 3.8 X 10-"^"* cm2 for a 70-GeV WIMP 
mass, the XENON-10 experiment has set a comparable 
limit of 8.8 X 10"'^'' cm^ for a 100-GeV WIMP mass [3, 
showing that liquefied noble gases are viable dark matter 
targets. 

Events in a noble liquid dark-matter detector may arise 
from scattering off of the nucleus or atomic electrons; 
dark matter will only scatter off the nucleus to an appre- 
ciable extent. The ratio of the scintillation light yield for 
nuclear recoil events relative to electronic recoil events is 
defined as the scintillation efficiency or Loff- 

A WIMP dark matter search requires an energy thresh- 
old on the order of tens of keV, and it is necessary to 
measure the scintillation efficiency down to this energy 
threshold so as to quantify the WIMP detection sensi- 
tivity. In order to make this measurement of Leff, a 
D-D neutron generator was used to produce neutrons 
that scattered from a liquid argon detector into an or- 
ganic liquid scintillator detector used as a coincidence 
trigger. The organic scintillator was placed at a series of 
known angles, and the energies of the selected nuclear re- 
coils in the liquid argon were kinematically determined. 
The scintillation efficiency was determined from the ra- 
tio of the measured electron-equivalent recoil energy at a 
given scattering angle to the expected nuclear recoil en- 
ergy (keVr) at that angle. Details of this measurement 
in a 4-kg liquid argon detector are presented in this pa- 
per, along with scintillation efficiency results for nuclear 
recoil energies between 10 and 250 keVr at zero electric 



field. 



II. REVIEW OF PHYSICAL PROCESSES AND 
MEASUREMENTS 

Discrimination between nuclear recoil events that char- 
acterize a WIMP signal and electronic recoil events 
that characterize the primary backgrounds is essential 
in WIMP detectors, particularly for the case of liquid 
argon which contains the radioactive isotope ^^Ar. The 
noble liquid detectors use two methods to achieve this 
discrimination. Single-phase detectors use pulse-shape 
discrimination (PSD) based solely on scintillation light 
to discriminate between event types, while dual-phase 
detectors can collect both scintillation light and ioniza- 
tion, employing a combination of PSD and the relative 
size of the light and ionization channels to identify events. 
PSD is made possible because ionizing radiation in liquid 
noble gases results in the formation of excited diatomic 
molecules (excimers) that can exist in either singlet or 
triplet states, with very different lifetimes. In liquid ar- 
gon these lifetimes are 7 ns and 1.5 /xs, respectively [olIToj 
and the scintillation light is produced in the decay of 
these states. As different types of excitation produce dif- 
ferent ratios of triplet to singlet molecules, the relative 
amplitudes of the fast and slow components can be used 
to determine what type of excitation occurred. The ef- 
fectiveness of this PSD is directly dependent on the num- 
ber of detected photoelectrons in an event, and thus the 
light yield for both nuclear recoils and electronic recoils 
sets the energy threshold for which electronic recoil back- 
grounds are negligible, in turn determining the ultimate 
sensitivity of the detector to dark-matter- induced nuclear 
recoils. 

The excimers that provide the scintillation light are 
formed in two ways. An excited atom (exciton) can com- 
bine with another atom in the liquid to produce the ex- 
cimer, or an ionized atom can combine with another atom 
to form a diatomic ion, which in turn recombines with 
an electron, eventually resulting in the production of the 
excimer. The ratio of exciton production to ion pair pro- 
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duction in liquid argon has been caleulated to be 0.21 
[Tl| . indicating that the majority of the scintillation light 
in liquid argon comes from excimers formed indirectly 
from argon ions, rather than directly from excited argon. 
The average energy required to produce an electron-ion 
pair in liquid argon has been measured to be 23.6 ± 0.3 
eV , and the average energy needed to produce a sin- 
gle photon has been calculated to be 19.5 ± 1.0 eV p^ . 
From this, the maximum possible scintillation yield in 
liquid argon is about 51 photons per keV of deposited 
energy, in the extreme case where the excimer formation 
and scintillation processes are perfectly efficient. 

In actuality, the absolute light yield is reduced through 
a number of different mechanisms. Energy may be lost 
by means other than exciton and ion pair formation, the 
excitons may undergo non-radiative collisions, and the 
recombination of diatomic ions may be incomplete. The 
first mechanism is known to be significant for nuclear re- 
coils, for which a significant portion of the e nerg y is lost 
to atomic motion as described by Lindhard |14l |. Thus, 
the scintillation light yield is expected to be reduced for 
nuclear recoil events compared to electronic recoils. Mea- 
surements of the scintillation efficiency for nuclear recoils 
relative to electronic recoils in liquid xenon [l5l - ll9| indi- 
cate that there is an additional reduction in the nuclear 
recoil scintillation yield due to collisions between free ex- 
citons that result in an ion and a ground-state atom, as 
described by Hitachi [2^ . The rate of these biexcitonic 
collisions is dependent on the density of the excitations, 
thus the amount of quenching increases with increasing 
linear-energy-transfer (LET) and is most significant at 
larger recoil energies. This mechanism of biexcitonic 
quenching is expected to apply to liquid argon as well 
as liquid xenon, and a model for scintillation efficiency 
in argon, neon and xenon taking LET into account has 
been proposed by Mei et al. [21| A further reduction in 
scintillation yield can result when some fraction of the 
ion-electron pairs do not recombine to produce an ex- 
cimer and the electrons escape instead [13|. 



III. EXPERIMENTAL DETAILS 

The liquid argon scintillation efficiency was measured 
using the MicroCLEAN detector at Yale University. The 
active volume is 3.14 liters of liquid argon viewed by 
two photomultiplier tubes (PMTs). Figure [T] shows a 
schematic of the central region and PMTs. The ac- 
tive region is defined by a teflon cylinder 200 mm in 
diameter and 100 mm in height, with two 3-mm-thick 
fused-silica windows closing the top and bottom. Two 
200-mm-diametcr Hamamatsu R5912-02MOD photomul- 
tiplier tubes (PMTs) arc held in place by tefion rings 
above and below the central volume and view the active 
region through the windows. Because liquid argon scin- 



tillates in the vacuum ultraviolet at 128 nm [2J], all inner 
surfaces of the teflon and fused silica are coated with a 
thin fllm of tetraphenyl butadiene (TPB) The TPB 
shifts the wavelength of the ultraviolet light to approxi- 
mately 440 nm so that it may pass through the windows 
and be detected by the PMTs. Both windows are coated 
with 0.20 ± 0.01 mg/cm^ of TPB and the teflon cyhndcr 
is coated with 0.30 ± 0.01 mg/cm^ of TPB. The teflon 
cylinder, windows and PMTs are all immersed directly 
in liquid argon, contained within a 25-cm-diameter by 
91-cm-tall stainless steel vessel. 
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Relative scintillation efficiencies in liquid argon have 
been measured for a number of different particle types. 
For heavy fission fragments with kinetic energy around 80 
MeV, the scintillation efficiency relative to 1-MeV elec- 
trons has been measured to be 0.21 ± 0.04 (2^. For al- 
pha particles, the scintillation efficiency has been mea- 
sured to be 0.9 for 6-MeV alphas relative to 1-MeV elec- 
trons [H] and 0.4 for 5.3-MeV alphas relative to 1.2-MeV 
electrons [1^. The observed triplet lifetime from the 5.3- 
McV alphas was 800 ns. which may indicate additional 
absorption due to impurities, but both values are plausi- 
ble given the expectation that the scintillation efficiency 
for alphas should fall somewhere between that of fission 
fragments and unity. A previous measurement of scintil- 
lation efficiency for nuclear recoils in liquid argon by the 
WARP collaboration gives 0.28 ± 10%, measured at 65 
keV average recoil energy . 



FIG. 1; Schematic representation of the scintillation volume 
and PMT orientation. 

The stainless steel vessel is housed inside a vacuum 
dewar, and liquid argon is introduced though a tube on 
the top of the vessel. The argon is liquefied from puri- 
fied gas in a copper vessel mounted to the end of a Cry- 
omech PT805 pulse-tube refrigerator. All components 
that come into contact with the gas or liquid are baked 
to at least 60° C, and the ultra-high-purity argon gas 
(99.999%) is passed through a heated Omni Nupure III 
gas-purification getter before entering the vessel. Out- 
gassing can cause impurities to build up in the detec- 
tor, decreasing the light yield by quenching the argon 
excimers or absorbing the UV scintillation photons. To 
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avoid signal degradation, the argon is continually circu- 
lated through the getter and reliquefied at a rate greater 
than 2.0 slpm. No reduction in signal was observed dur- 
ing the run. PMT signals were ditigized using an 8-bit 
500 MSPS waveform digitizer with each of the PMTs 
capturing both low-gain and high-gain waveforms. More 
details about the experimental apparatus, data acquisi- 
tion, and purity measurements are available in fTo|. 
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FIG. 2: Example of an electronic recoil-induced scintillation 
event in liquid argon. 

A sample oscilloscope trace from an electronic recoil 
scintillation event in argon is shown in Fig. [5] A 10- 
/iCi sealed ^''Co source is used for daily measurements of 
the scintillation light yield for electronic recoils, with a 
sample spectrum shown in Fig. [3] along with results from 
the simulation to be discussed in the next section. This 
source produces 122-keV, 137-keV and 14.4-keV gamma- 
rays with branching ratios of 86%, 11% and 9%, respec- 
tively. Spectra were taken for each day of data taking, 
with a Gaussian fit to the 122-keV peak providing a scin- 
tillation signal yield calibration for that day in units of 
photoelectrons per keV of energy deposited by an elec- 
tronic recoil, denoted as photoelectrons/keVee. Over the 
course of the four-month run, the signal yield remained 
stable to within 5% at 4.85 ±0.01 photoclectrons/keVee. 
To check the quality of the energy calibration, a 10-^Ci 
^^Na source that produces 511-keV gamma-rays is used 
as a second point of reference, and the 511-keV line ap- 
pears at a photoelectron yield that is within 1% of the 
value predicted from ^^Co source calibration. 

A portable Thermo Electron MP320 deuterium- 
deuterium (D-D) neutron generator is used as a neu- 
tron source, with an organic scintillator detector as a 
secondary coincidence trigger. The experimental setup 
can be seen schematically in Fig. |4l In the forward di- 
rection, the D-D generator produces 2.8-MeV neutrons. 
Some of these neutrons scatter in the liquid argon, and 
for a given position of the organic scintillator, only neu- 
trons that scatter at a specific angle are selected by the 
coincidence trigger. By changing the angle the organic 
scintillator makes with the neutron generator-liquid ar- 
gon detector axis, the scattering energy of the recoil nu- 
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FIG. 3: (Color online) Plot of the *^Co spectrum, along with 
a simulation of the expected spectrum. 



cleus in the liquid argon can be varied according to the 
following equation: 



COS 



[m„ 4- niAr — "n^n cos^(0) — 



(1) 



where Ein is the incident neutron energy (2.8 MeV), A 
is the atomic mass number, and 9 is the scattering angle 
of the outgoing neutron. Data were taken at 19 different 
angles corresponding to recoil energies between 10 keV 
and 250 keV. 
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FIG. 4: (Color online) Top view of the neutron scattering 
setup. Shown are the neutron generator and the organic scin- 
tillator. The size of the argon cell is not representative. 



IV. MONTE CARLO SIMULATIONS 

In order to understand the data, we developed a Monte 
Carlo simulation of the argon detector, cryostat, organic 
scintillator and surrounding lab space. The software 
framework used was RAT, which combines Geant4f26j. 
CLHEPfl^l, and ROOTjlil into a single simulation and 
analysis package. A detailed optical model of the in- 



ner detector and PMTs is included in the Monte Carlo 
which allows us to estimate smearing of the detected sig- 
nal. While this model gives results that are in fairly good 
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agreement with our gamma calibrations, we add an ad- 
dition smearing term for the neutron scattering analysis 
to take into account the lower photon yield for a given 
energy. 

An image of the detector geometry can be seen in 
Fig.m The argon detector is in the central vertical cylin- 
der, where the various layers of steel are set to be semi- 
transparent so that the inner workings are visible. While 
there was only one organic scintillator detector at a given 
time in the real experiment, the simulation used many 
such detectors, so as to allow Monte Carlo data for half 
of the organic scintillator positions to be collected simul- 
taneously and reduce computation time. In the picture, 
the cylinders off to the right of the argon detector repre- 
sent the organic scintillator in its various positions. As 
the adjacent positions of the organic scintillator overlap, 
independent simulations were performed for each set of 
colored organic scintillator positions. 




FIG. 5: (Color online) Monte Carlo detector geometry. 



In addition to the neutron scattering simulation, the 
detector response to an external '^^Co gamma source was 
also modeled. The origin of the 122 keV and 137 keV 
gamma rays was set just outside of the outer vacuum 
can as in the real detector. The results of this simulation 
can be seen in Fig. [3] showing very good agreement with 
the experiment. 



V. ANALYSIS 

The analysis begins by combining the high-gain and 
low-gain waveforms for each PMT in a given event. To 
do this the first two microseconds of each PMT's high- 
gain and low-gain waveforms are separately averaged to 
calculate baselines. These baselines are then subtracted 



to give zero-offset waveforms. Each high gain waveform 
is then scanned to determine if the waveform digitizer is 
saturated. In the case of saturation, the high-gain and 
low-gain waveforms are aligned and the low-gain wave- 
form's samples are inserted where the high-gain waveform 
was saturated. 

At this point each PMT waveform is scanned to de- 
termine the trigger time of the event. The trigger time 
for each PMT is defined as the time where the wave- 
form reaches 20% of its maximum height. The average 
of the two trigger times is taken as the start time for the 
event and a timing cut is applied to the two PMT wave- 
forms to remove events where the difference in timing is 
greater than 20 ns. The waveforms from both PMTs are 
then integrated in two timing intervals, the first from 20 
ns before the trigger time to 100 ns after and the sec- 
ond from 100 ns to 5 fis. The region between 5 fis and 
lA fis is scanned for single photoelectron pulses and used 
to determine the single photoelectron spectrum. Any 
region where the waveform's voltage value exceeded ap- 
proximately one third of a single photoelectron's peak 
voltage is integrated from 10 ns preceding the crossover 
sample to 50 ns following it. After this procedure has 
been performed on every event in a run, the run's single 
photoelectron value is fit and used to convert the inte- 
grated waveform charges into photoelectrons. 

A PMT asymmetry cut is used to remove events that 
are near the windows of the detector. The asymmetry 
is defined as the difference in the signals observed by 
the two PMTs divided by their sum, and events with an 
asymmetry of more than 60% are removed. Events with 
approximately 2000 photoelectrons or more can cause the 
saturation of one or both of the detector PMTs, and since 
this will cause events to have poor energy reconstruction, 
a cut is applied to remove events in which either PMT's 
output becomes greater than 2 V. This cut removed nu- 
clear recoil events with energies above 110 keVee for runs 
with recoil angle below 125° degrees and nuclear recoil 
events with energies above 180 keVee for 125° and 142° 
degree runs. All cuts applied to the data up to this point 
are considered to be data quality cuts and were applied 
to both ^^Co and neutron runs. 

To distinguish neutron scatters from other back- 
grounds in the neutron data sets, two additional cuts 
involving the organic scintillator are applied. The first 
is a time-of-flight (TOF) cut which removes events in 
which the organic scintillator is triggered before the de- 
tector, as well as events in which the neutron arrives late 
due to multiple scatters. In addition, this cut also helps 
remove background gammas from our neutron data sets. 
The position and width of this cut is set by the location 
of the single scattering neutrons in the Monte Carlo TOF 
spectrum. 

The second cut associated with the organic scintillator 
uses the pulse shape in the organic scintillator to distin- 
guish between neutrons and background events. We use a 
scatter plot of the pulse area within 100 ns of the organic 
scintillator trigger versus the organic scintillator wave- 
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form's maximum voltage, shown in Fig[6l to determine a 
quadratic curve that separates two types of events. This 
curve divides the scatter plot into two distinct regions: 
an electronic recoil band that appears for all types of runs 
and a nuclear recoil band that is only present in neutron 
runs. The separation of these regions did not appear to 
change with recoil energy and variation of the regions 
had a negligable effect on the results. 




Charge (pC) 

FIG. 6: (Color online) An example of the pulse shape cut 
applied to the organic scintillator. A scatter plot of the pulse 
area versus peak pulse height shows a distinct region of neu- 
tron events bounded by the red quadratic curve and lines. 



We apply one final cut to remove electronic recoils 
from the sample, exploiting the PSD to discriminate be- 
tween event types in liquid argon. Based on our previous 
work [l^l, we define a discrimination parameter, Fprompt, 
as the fraction of light arriving in a prompt time window. 
For comparison, the mean Fprompt for recoils between 
5 and 32 kcVce range from 0.39 to 0.28 for electronic re- 
coils and 0.56 and 0.7 for nuclear recoils [l^l. We apply a 
relatively loose cut, removing events with Fprompt < 0.35 
from the final sample. We checked the analysis with a 
tighter cut of Fprompt < 0.50 and incorporated the dif- 
ferences in the systematic uncertainties. 

To determine the scintillation efficiency of the liquid 
argon, the measured energy spectrum for each scatter- 
ing angle is compared to the energy spectrum produced 
by the Monte Carlo simulation analyzed with the same 
asymmetry and TOF cuts. We also use the Monte Carlo 
simulation to calibrate the TOF cut. For each scattering 
angle in the Monte Carlo simulation, a TOF spectrum 
is produced using events that only scattered once in the 
detector before reaching the organic scintillator. This al- 
lows us to find the range of the TOF for single scattering 
neutrons for each position. After applying the same TOF 
and PMT asymmetry cut to the simulation as used for 



the data, a Monte Carlo energy spectrum is then gener- 
ated for each recoil angle. 

There are two convolutions applied to the Monte Carlo 
recoil spectra before fitting them to the data. First, 
to account for the variation in the single photoelectron 
charge, the Monte Carlo photon counts are smeared 
using the measured single photoelectron charge distri- 
bution from the photoelectron calibration data. Sec- 
ond, since the simulations were performed assuming a 
100% scintillation efficiency, an additional smearing of 
3.25 X -^(I — Lctf)Npc, is applied to the Monte Carlo to 
account for the difference in counting statistics between 
a scintillation efficiency of 100% and that obtained from 
the data. The proportionality constant of 3.25 empir- 
ically account for the observed broadness of the clearly 
resolved peaks at 191, 211 and 239 keVr scattering angles. 
It is well known that noble liquid detectors do not reach 
the ideal energy resolution predicted by Poisson photo- 
electron statistics, largely due to ionization-scintillation 
anti-correlation[l|, [2^ . 

The MINUIT fitting package is used to perform a fit 
of the Monte Carlo to the data with the normalization 
and the scintillation efficiency as free variables. Each 
Monte Carlo spectrum is binned using the same binning 
as the corresponding recoil data and used to generate a 
spline for fitting. First, the entire range of the data is 
used in the fit. Then, a Gaussian is fit to all events in 
the Monte Carlo identified as singly-scattered neutrons. 
This fit is used to define a new fit range consisting of 
plus and minus three sigma around the centroid of the 
Monte Carlo single scattered neutron distribution. The 
final fits are performed over a restricted range around the 
single scattering Monte Carlo neutron distribution where 
we expect to observe our signal. This fitting procedure 
was checked using the Monte Carlo sample with a set 
Lcs value of 0.25 and was able to recover this set Leff at 
each recoil angle. 

The results for all organic scintillator positions are pre- 
sented in Fig. [7l and the scintillation efficiency as a func- 
tion of energy is shown in Fig. \8\ After studying the sys- 
tematic effects described in the next section, we found 
that the individual scintillation efficiency values were 
constant across the range of recoil energies studied above 
20 keVr, with a mean of 0.25 ± 0.02 4- O.Ol(corrclated). 
Below 20 keVr, our data exhibit an upturn in scintilla- 
tion efficiency as the energy goes to zero, and we were 
unable to find an experimental cause for this upturn. It 
is therefore unknown if this is a physically real effect or 
if we lose our ability to distinguish nuclear recoils from 
other backgrounds at these low energies. All observed 
values and uncertainties are listed in Table HI 



VI. SYSTEMATIC UNCERTAINTIES 

The sources of uncertainty considered are categorized 
as those associated with detector operation, triggering ef- 
fects, Monte Carlo background normalization, TOF win- 
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FIG. 7: (Color online) Plotted are the recoil energy spectra for the 19 organic scintillator positions used in this experiment. 
The data are taken with the organic scintillator located at the angle indicated in the legend of each plot with the corresponding 
recoil energy indicated just below. The fit value for Leff with the statistical uncertainty from the fit is also listed in the legend 
of each plot. In each plot, the upper (red) histogram is the output of the GEANT4 based Monte Carlo simulation of single and 
multiple neutron scatters in the detector. This upper histogram is fit to the data in the solid region, whereas the dotted part 
shows the MC simulation outside of the fit range. The lower solid histogram (blue) is the subset of the Monte Carlo events 
where the neutron only scatters once in the detector volume. 
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FIG. 8: (Color online) Scintillation efficiency as a function of 
energy from 10 to 250 keVr. The weighted mean (red line) 
is generated from the data above 20 keVr and puts the mean 
scintillation efficiency at 0.25. The value measured by WARP 
is 0.28 at 65 keVr 0]. 



dow and fit range effects. The uncertainties from the 
sources discussed in this section are all combined and 
the final resulting uncertainty for each scattering angle 
can be found in Table. HI 

The first group of considered uncertainties deals 
with the data taking and stability of the neutron and 
^^Co data runs. Since the ^"^Co runs are used to cali- 
brate the light yield of the detector, the fit error of the 
^''Co peak and its stability over time directly affect the 
measured scintillation efficiency. These are estimated to 
be 2% and 1.6% respectively. There is a second uncer- 
tainty in the angle of the organic scintillator relative to 
the neutron generator that in turn introduces an uncer- 
tainty in the corresponding energy via Eq. [TJ We have 
determined the uncertainty of the angular position of the 
organic scintillator to be 1.3° at each position. 

We examined the effects of the trigger efficiency, specif- 
ically looking to address the upturn observed at low ener- 
gies which could be explained by a bias introduced by the 
trigger level. We took data for the 22° run at three dif- 
ferent hardware triggers and the 26° run at two different 
hardware triggers, and we examined the effect of hard- 
ware and software triggers on both the asymmetry cut 
and the final scintillation efficiency values. In all cases, 
the scintillation efficiency distributions did not system- 
atically change by varying the cuts and hardware thresh- 
old. We also performed a toy Monte Carlo using the 
time dependence of the scintillation light [TO] and the ob- 
served single photoelectron distribution to estimate pos- 
sible threshold effects. This study found the effect of 
any threshold bias given our hardware trigger level to be 
less than 1%, much smaller than the other errors in the 
measurement. Therefore, we conclude that the triggering 
threshold does not explain the upturn at low energies. 

A third source of uncertainty arises because the Monte 
Carlo simulation does not exactly reproduce the observed 



background shape, as is evident from comparing the 
dashed and solid red lines in the large angle scattering 
spectra of Fig. [T] To account for this inconsistency, the 
data for each recoil energy were reanalyzed under the 
assumption that the size of the multiple scattering back- 
ground in the histograms used to perform the fits varied 
by ±50% relative to that predicted by the Monte Carlo. 
The variations observed in this reanalysis arc approxi- 
mately 10% below 20 keVr and 2% above 20 kcVr and 
are included in the errors listed in Table HI 

To determine the uncertainty due to the TOF cuts, the 
TOF window was expanded separately up and down in 
time by 50%. This allows for recoil neutrons with smaller 
TOFs to be included when the window is expanded down- 
ward and larger TOF neutrons when expanded upward. 
The effect of this variation was mostly in the lowest three 
data points and allowed them to move downwards in scin- 
tillation efficiency by about 0.04. 

There is an uncertainty associated with fitting the data 
in a limited range around the predicted single-scattered 
neutron peak position. To estimate this uncertainty, we 
expand the fit range to include ±5 sigma around the ccn- 
troid of the single-scattered neutrons, instead of the 3- 
sigma range used in the standard fit. The result of the 
wider fit range is to systematically push the determined 
scintillation efficiency up for energies between 20 keVr 
and 120 keVr. This effect appears to be caused by a 
disagreement in the high energy tails of the data and 
Monte Carlo, similar to the disagreement observed be- 
tween Monte Carlo and data at lower energies for the 
high recoil angles. Changing the fit range adds a corre- 
lated error of -1-0.01 to the measured mean scintillation 
efficiency in Fig. IS] 



VII. RESULTS 

The scintillation efficiency for nuclear recoils in liq- 
uid argon has been measured relative to electronic re- 
coils for nuclear recoil energies from 10 keVr to 250 
keVr. The scintillation efficiency values found at each 
of the recoil angles can be found in Table |T] and are 
plotted in Fig. |51 The ratio of the nuclear recoil scin- 
tillation response to the electronic recoil response is 
0.25 ± 0.02 + 0.01 (correlated) for recoils above 20 keVr. 
The bin to bin correlation coefficient of the scintillation 
efficiency for the nuclear recoils above 20 keVr was 0.41. 
An observed upturn in the scintillation efficiency below 
20 kcVr is currently unexplained. The scintillation effi- 
ciency for nuclear recoils should also have some depen- 
dence on any applied electric fields and this is a topic for 
further research. 
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TABLE I: Table of energies and scintillation efficiencies from 
Fig. |8] LcH values for energies above 20 keVr also have an 
additional correlated error of +0.01. The uncertainties shown 
are the combined statistical and systematic uncertainties. 
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